function obj=obj_y(y, rho, kappa, zeta, logit_coef, X)


y=exp(y); % make sure y is positive
obj=-exp_u(y, rho, kappa, zeta, logit_coef, X);

end